Accepted to ApJ, August 6 2009; published October 1 2009 

Preprint typeset using LSTjhX style emulateapj 



DEEP KECK U -BAND IMAGING OF THE HUBBLE ULTRA DEEP FIELD: 
A CATALOG OF Z - 3 LYMAN BREAK GALAXIES 



Marc Rafelski, 1 Arthur M. Wolfe, 1 Jeff Cooke, 



Hsiao-Wen Chen, 3 Taft E. Armandroff, 4 & Gregory D. Wirth 4 



O 
O 

(N 
CU 

(N 

o 

u 

43 

Oh- 
i < 

o ■ 



> 

m 
q 

06 
o 

On 
O 



13 



Accepted to ApJ, August 6 2009; published October 1 2009 

ABSTRACT 

We present a sample of 407 z ~ 3 Lyman break galaxies (LBGs) to a limiting isophotal M-band magnitude of 
27.6 mag in the Hubble Ultra Deep Field (UDF). The LBGs are selected using a combination of photometric 
redshifts and the M-band drop-out technique enabled by the introduction of an extremely deep M-band image 
obtained with the Keck I telescope and the blue channel of the Low Resolution Imaging Spectrometer. The 
Keck M-band image, totaling 9 hrs of integration time, has a la depth of 30.7 mag arcsec -2 , making it one of 
the most sensitive M-band images ever obtained. The M-band image also substantially improves the accuracy 
of photometric redshift measurements of ~ 50% of the z ~ 3 LBGs, significantly reducing the traditional 
degeneracy of colors between z ~ 3 and z ~ 0.2 galaxies. This sample provides the most sensitive, high- 
resolution multi-filter imaging of reliably identified z ~ 3 LBGs for morphological studies of galaxy formation 
and evolution and the star formation efficiency of gas at high redshift. 

Subject headings: cosmology: observations — galaxies: distances and redshifts — galaxies: evolution — 
galaxies: general — galaxies: high-redshift — galaxies: photometry 



1. INTRODUCTION 



The Hubble Ultra Deep Field (UDF: iBeckwith et al.| [2006) 
provides the most sensitive high-resolution images ever taken, 
yielding a unique data set for studying galaxy evolution. 
These data have contributed to many scientific advances, in- 
cludin g constraining the st ar formation efficiency of gas at 
z ~ 3 (IWolfe & Chen 2006), aiding in determ ining the lumi- 
nosity function in the redshift range 4 < z < 6 dBouwens et alJ 
2007), bringing insig ht into the merger fractions of galaxies 
dConselice et al.ll2008l) . a nd yielding the discove ry of clumpy 
galaxies at high redshift (Elmegreen et al. 2003). 

Knowledge of galaxy redshifts is essential to understand- 
ing their nature, and various approaches are used to esti- 
mate this key attribute. A number of studies identify objects 
in the redshift range 4 < z < 6 using the deep multi-filter 
(BVIZJH) UDF images by identifying so-called "dropout" 
galaxies which are detectable in certain broadband filters 
but no t in others (IBeckwith et al.ll2006t iBouwens et~aUl2006L 
2007). Others use photometric redshifts derived across the 
entire redshift range < z < 6 by analyzing the co lors of 
galaxies in a wider range of filters dCoe et a l. 2006). This 
latter approach has t he potential to provide rela tively accu- 
rate (cfaz/(i+z) ^5 0-1, Fernan dez-Soto et aTJ l2~001) redshift es- 
timates for a large number of galaxies in a given field, but 
traditionally has serious problems producing accurate results 
near z ~ 3, a redshift range for which key spectral energy dis- 
tribution (SED) features fall blueward of the previously avail- 
able filter set. 

The majority of galaxies observed at z ~ 3 are Lyman break 
galaxies (LBGs): star-forming galaxies that are selected based 
on the break in their SED at the 912 A Lyman limit primar- 
ily by interstellar gas intrinsic to the galaxy, as well as a flux 
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decrement shortward of 1 2 1 6 A in the galaxy rest frame due to 
absorption by the Lyman series of opticallythick hydrogen gas 
along the line of sight. This Lyman limit discontinuity in the 
SED significantly dims these galaxies shortward of ~ 3500 A, 
allowing them to be found with the U-band dropout te chnique 
dSteidel & Harmltori[l992T: ISteidel et al J[T995l 1 1 996alfbh . Pre- 
viously, the available observations in the UDF did not include 
deep M-band imaging. The next prominent broadband signa- 
ture in the SED is the 4000 A break, which is redshifted to 
the near infrared (IR) for z ~ 3 LBGs. Without the M-band or 
very deep IR coverage, it is very difficult to determine from 
broadband imaging alone whether the observed decrement in 
the SED near the observed-frame ~ 4800 A is a low-redshift 
galaxy with a 4000 A break, or a high-redshift galaxy with a 
decrement from the 1216 A break. This degeneracy causes 
"catastrophic" errors in the photometric redshifts of galaxies 
atz ~ 3 without M-ba nd data dEllisll9 97: Fernande z-Soto et"aTJ 
fT^lB^nfte3l2000h . 

The purpose of this paper is to present a reliable sample of 
LBGs at z ~ 3 in the UDF through the introduction of ultra- 
deep u'-band imaging acquired with the 10m Keck I tele- 
scope. The Keck I tel escope and the Low Resolution Imagin 
Spectrometer (LRIS; lOke et al][l99l iMcCarthv et aTJ[T99^ 
form an ideal combination to probe galaxies at the z ~ 3 epoch 
due to the light-gathering power of the 10 m primary mir- 
ror and the outstanding efficiency of the LRIS blue channel 
in the near UV. This allows the M-band filter (with effective 
wavelength A - 3400 A, FWHM- 690 A) used with the 
blue arm (LRIS-B) to be ~ 300 A bluer and ~ 360 A wider 
than the M-band filter on the Visible Multi-Object Spectro- 
graph (VIMOS: iFevre et al.ll2003l) instrument (A c - 3700 A, 
FWHM~ 330 A) at the Very Large Telescope (VLT) and still 
be effective. Consequently, this enables us to probe the Ly- 
man break efficiently to a lower limit of z ~ 2.5, versus the 
VIMOS limit of z ~ 2.9 (a s shown by the VIM OS observa- 
tions of GOODS-South by iNoninoetai] 120091) . In addition, 
although the UDF field can only be observed at high air mass 
from Mauna Kea, the total throughput of LRIS-B and its bluer 
M-band is approximately twice that of VIMOS with its redder 
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M-band. This gives Keck the unique ability to select lower 
redshift LBGs via their Lyman break. 

We assemble a reliable sample of LBGs in the UDF with a 
combination of the M-band dropout technique and photomet- 
ric redshifts, and provide photometric redshifts for all galaxies 
with good w-band photometry. We find that the M-band imag- 
ing improves the photometric redshifts of z ~ 3 galaxies by 
reducing the degeneracy between low and high-redshift galax- 
ies. Furthermore, the combination of deep LRIS M-band and 
high-resolution multiband Hubble Space Telescope (HST) 
imaging provides an unprecedented view of z ~ 3 LBGs to 
improve our underst anding of their h ighly irregular rest-frame 
UV morphologies (lLaw et al]l2007l) down to unprecedented 
depths. 

The deep, high-resolution LBG sample will also extend 
current constraints on the star formation efficiency of gas at 
z ~ 3. Reservoirs of neutral gas are needed to provide the 
fuel for star formation. Damped Lya systems (DLAs), se- 
lected for their neutral hydrogen column densities of Nm > 
2 x 10 20 cm~ 2 , dominate the neutral-gas content of the universe 
in the redshift interval < z < 5. DLAs contain enough gas to 
account for 50% of the mass conten t of visible matter in mod- 
ern galaxies (see lWolfe et alJ l2005. for a review) and may act 
as neutral-gas reservoirs for star formation, since stars form 
when local values of the H I column density, Nm, exceed a 
critical value. At high redshift, the star formation rate (SFR) 
is assumed to follow the Kennicutt-Schmidt (KS) law (estab- 
lished for nearby galaxies) which states how the SFR per unit 
physical area, E*, relates to the neutral gas (i.e., H I and H2) 
column density: E* oc N^ 4 (lKennicuttlll998h . Strong DLAs 
with N m > 1 .6 x 10 21 cm" 2 in the redshift range 2.5 < z < 3.5 
are predicted to have emission from star formation in the rest- 
frame FUV redshifted into the optical such that they are de- 
tectable in the UDF F606W image. IWolfe & Chenl d2006h 
searched for low surface-brightness emission from DLAs in 
the UDF and found the in situ SFR efficiency of DLAs to 
be less than 5% of the KS law. In other words, star forma- 
tion must occur at much lower rat es in DLAs at z ~ 3 than 
in modern galaxies. Whereas the Wolfe & Chenl d2006l) re- 
sults set sensitive upper limits on in situ star formation in 
DLAs excluding known galaxy regions, no such limits exist 
for DLAs containing LBGs. The sample presented here en- 
ables a search for spatially extended low-surface-brightness 
emission around z ~ 3 LBGs which will yield constraints on 
the star formation efficiency at high redshift (Rafelski et al. 
in prep). This is one of the main motivations for constructing 
this sample, and therefore we construct our sample conserva- 
tively to minimize the number of potential interlopers. 

We present our catalog of z ~ 3 LBGs, provide photomet- 
ric redshifts for the entire sample of objects that have reliable 
M-band photometry, and make the M-band image available to 
the public. The observations are described in §2, and the data 
reduction and analysis in §3. We discuss the photometric se- 
lection of z ~ 3 galaxies in §4, and summarize our major find- 
ings in §5. Throughout this paper, we adopt the AB mag- 
nitude system and an (H,M,^A,h) = (0.3,0.7,0.7) cosmology 
with parameters T = 0.21,n = 1, which are largely consistent 
with recent values (Hin shaw et al.ll2009h . 

2. OBSERVATIONS 

The M-band (A c - 3400 A, FWHM- 690 A) images of 
the UDF (a(/2000) = 03 h 32 m W, (5(72000) = -27°47'29."1) 
were obtained with the 10 m Keck I telescope using the LRIS. 
The M-band data were taken with the new Cassegrain Atmo- 



spheric Dispersion Corrector (Cass ADC; Phillips et a U2006l) 
to minimize image distortions from differential atmospheric 
refraction. The Cass ADC was critical to the success of these 
observations because of the low elevation of the UDF from 
Mauna Kea, and the blue wavelength of our primary band. 
The data also benefit from the backside illuminated, dual UV- 
optimized Marconi 2048 x 4096 pixel CCDs on LRIS-B, with 
0."135 pixels and very high UV quantum efficiency (~ 50% 
at A ~ 3450 A). Because the quantum efficiency of the two 
CCD chips varies ^30%-35% in the M-band, we placed the 
UDF entirely on the more sensitive chip (CCD1). We used a 
dichroic beam splitter (D460) to simultaneously observe the 
M-band on the blue side, and the V-band and /?-band on the red 
side. The red channel data were taken for astrometric and pho- 
tometric consistency checks, and were not used in this study 
other than for calibration purposes due to the much deeper 
and higher resolution UDF images available over those wave- 
length ranges. 

The observations were carried out in darktime over two 
runs, 2007 October 7-9 (three half nights) and 2007 Decem- 
ber 3-4- (two half nights). We lost the entire first night of 
the first run to weather, and had moderate weather and seeing 
conditions throughout both runs that yielded a median seeing 
FWHM of ~1."3 in the M-band. In ord er to maximize time on 
sky, w e adopted the dither strategy of Sawi cki & Thom pson 
(2005, see their Eqn. 1), setting the red channel exposure 
times such that the last readout of the red channel coincided 
with the end of the blue channel readout. The M-band images 
were acquired as a series of 36 x 900s exposures to avoid the 
nonlinear regime of the CCD, totaling 9 hrs of integration time 
on target. We executed a nine-point dither pattern with 10" 
dithers to deal with bad pixels and to create a super-sky flat. 

The UDF can only be observed at large air masses at Mauna 
Kea, which can affect the shape of the M-band throughput. We 
show a histogram of our observed air masses and the variation 
of the filter throughput for different air masses in Figure [T] 
We find that the variations in air masses in our sample do not 
significantly affect the blue side cutoff of our filter. Specif- 
ically, we find that the variation between our best and worst 
air mass yields a change in X„ < 10, with typical changes in 
A e being < 5. Given the small variation of A G , we derive the 
final M-band filter transmission curve by convolving the mea- 
sured filter throughput with the atmospheric attenuation at our 
average air mass of 1 .57 and the CCD quantum efficiency. 

Throughout the paper we utilize the B, V, i', and z' band 
(F435W, F606W, F 775W, and F850LP, r espectively) obser- 
vations of the UDF (Bec kwith et alJl2006h . obtained with the 
Wide Field Ca mera (WFC) on th e HST Advanced Camera for 
Surveys (ACS: iFord et alJl2002h . These images cover 12.80 
arcmin 2 , although we prune our catalog to the central 1 1 .56 
arcmin 2 which contains at least half the average depth of the 
whole image and overlaps our M-band image with uniform 
depth. In addition to these ACS images, we also include ob- 
servations taken with the NIC MOS camera NIC3 in th e J and 
H bands (Fl 10W and F160W: lthompson et alj|2006h . These 
red wavelengths cover the central 5.76 arcmin 2 of the UDF, so 
we only use them whenever the field of view (FOV) overlaps. 
Figure [2] plots the total throughput of the filters used in this 
paper: the Keck LRIS-B M-band, the HST ACS B, V, i', and 
z' bands, and the HST NICMOS J and H bands, and include 
the CCD quantum efficiency and atmospheric attenuation. 
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FIG. 1 . — Transmission of the u-band filter used in this study for different air masses (X). The filters are corrected for the atmosphere attenuation at different 
air masses and the CCD quantum efficiency of LRIS-B on the Keck telescope. The histogram shows the range of air masses in the study, with the minimum and 
maximum air mass being plotted on the transmission curve. The change in the effective wavelength A„ between air masses is typically < 5 A. 




FIG. 2. — Transmissions of the filters used in this paper, corrected for 
the CCD quantum efficiency, and in the case of the w-band, the atmosphere 
attenuation at the average air mass of 1.57. The ii-band is from LRIS-B on 
the Keck telescope, the B, V, i', and z' bands from the WFC on the HST ACS, 
and the / and H bands from the NICMOS camera on HST. 



3. DATA REDUCTION AND ANALYSIS 

3.1. Image processing 

The LRIS-B data were processed in a combination of cus- 
tom code in IDL, and standard data reduction algorithms from 
IRAF 5 . The images were bias subtracted, first from the over- 
scan region and then from separate bias frames to remove any 
residuals, before being trimmed to remove any vignetted re- 
gions of LRIS. Super-sky fiats were created using IRAF from 
the median of all the unregistered images with sigma clip- 
ping to remove objects and cosmic rays, which were then 
used to flat-field the images. These prove superior to dome 
and twilight sky flats in determining the CCD response in 
the w-band because of the short wavelengths, and yield ex- 
cellent flats. Dithering off sets were determined using cus- 
tom code and SExtractor dBertin & Arnoutsl 19961) to locate 

5 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 



bright objects common to all the images. In order to drizzle 
dFruchter & Hookl l2002) only once and keep correlated noise 
to a minimum, the images were distortion corrected using the 
solution provided by J. Cohen & W. Huang (private commu- 
nication, February 2008) and shifted to correct for the dither- 
ing offsets all at once using the geotran package in IRAF. 
We drizzled with pixf rac = 0.5 to improve the point-spread 
function (PSF) and set the pixel scale such that the pixels are 
integer multiples of the UDF pixels, 0."12, for reasons ex- 
plained in §3.4. To maximize the signal to noise (S/N), the 
images were weighted by their inverse variances. The drizzled 
images were then combined using the IRAF task imcombine, 
with sigma clipping to remove bad pixels and cosmic rays. 
These images were trimmed to include regions of uniform 
depth (11.56 arcmin 2 ), normalized to an effective exposure 
time of Is, and then background subtracted using the global 
background determined with SExtractor. An astrometric 
solution was applied with the IRAF task ccmap by matching 
bright and nearly unresolved objects in the UDF to those in 
the final stacked image. The final rms astrometric errors are 
between 0."02 and 0."03, negligible in comparison to the 1 ."3 
FWHMoftheM-bandPSF. 

3.2. Photometric Calibration 

Moderate weather yielded no completely photometric night 
for our calibration, and we therefore calibrate d to the Multi- 
wave length Survey by Yale-Chile (MUSYC; iGawiser et alJ 
2006), which covers the same part of the sky as our primary 
observations. We include all our filters (u, V, and R) for 
th is calibration. All our calibrations are in the AB95 system 
of iFukugita et alJ (119961) . hereafter referred to as AB magni- 
tudes. We used the IRAF tasks phot and f itparams to solve 
for zero-point magnitudes, air-mass correction coefficients, 
and appropriate color correction coefficients. Specifically, we 
use the equation: 

m = -2.51og 10 (F) + Z-cX-Y, (1) 

where F is the flux in counts/s, Z is the zero-point magni- 
tude, c is the air-mass coefficient, X is the air-mass, and Y 
is the color term. We allow a color term to account for any 
differences between the f/-band filter used by MUSYC and 
the M-ba nd filter in our observ ations, similar to the color term 
used by IGawiser et alJ d2006l) to correct for their differences 
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compared to the Johnson-Cousins filter set. The galactic ex- 
tinction of 0.0384 is subtracted from the zero-point magni- 
tude, using the relat ion A(u) = 4.8£'(B- V) i nterpreted from 
Carde lli et all d 19891) . where E(B-V) = 0.008 dBeckwith et all 
2006). The final results for the M-band calibration are a zero- 
point magnitude Z = 27 .80 ± 0.03, an air-mass coefficient term 
c = 0.41, an average air-mass X = 1.57, and a color term 
Y = (0.13 ±0.02) x (U-B) AB . We double check our cal- 
ibration using multiple observed photometric standard stars 
(lLandoltlll 992) over a range of air masses, and the zero-point 
and air-mass correction coefficients are consistent with those 
found when calibrating to the MUSYC catalog. As a result, 
we are confident that the M-band image is well calibrated. 

3.3. Depth of the u-band Image 

It is useful to characterize the depth of the M-band image, 
however, different definitions exist to describe the sensitivity 
of an image. Two commonly quoted limits are presented here: 
a measurement of the sky fluctuations of the image, and a 
limiting magnitude corresponding to a 50% decrease in object 
counts through Monte Carlo simulations. 

The sky noise of the image is measured via the pixel to pixel 
rms fluctuations in the image, best measured by fitting a Gaus- 
sian to the histogram of all pixels without sources. Sources are 
identified with the program SExtractor, with the threshold 
set such that the negative image has no detections. This yields 
a depth of 31.0 mag arcsec" 2 , \a u sky fluctuations. How- 
ever, since the image is drizzled, correlated noise between 
the pixels is introduced. The theoretical inc rease in noise due 
to pix f rac = 0.5 using equation 10 from iFruchter & Hookl 
(2002) is 20%. Alternatively, the correl ated noise can be esti- 
mated empirically with equation 2 from Fernande z-Soto et"ai1 
( 1999), which uses a covariance matrix to determine a small 
overestimate of the real error 6 . This results in a slightly more 
conservative depth of 30.7 mag arcsec -2 , la„ sky fluctuations, 
which is what we quote here. 

Additionally, to get a better sense of the usable d epth of the 
image, a limiting magnitude is ofte n quoted (e.g.. lChen et al.1 
l2002t ISawicki & Thompson] 12005). We define u\ im as the 
magnitude limit at which more than 50% of the objects are 
detected. The best way to determine u\i m is through Monte 
Carlo simulations, which take into account both the sky sur- 
face brightness and the seeing in our image. Since our image 
PSF has more flux in the wings than a Gaussian, we plant 
both Gaussian objects, and objects modeled to fit our PSF 
using a two-dimensional (2-D) Moffat profile 7 . The custom 
IDL code was used to extract bright unresolved objects in 
the M-band image, take the median of all these objects, and 
create a composite object stack. The composite image was 
then fit bot h by a Gaussian, a nd by a 2-D Moffat profile us- 
ing MPFIT (Markwardt 2009), with a modification to ensure 
that the wings of the profile go to for the Moffat profile. 
We semi-randomly insert these objects with a range of fluxes 
into the M-band image. The locations of the planted objects 
are constrained such that they do not: 1) fall off the edges, 
2) fall on a real detected object, and 3) fall on any previously 
planted objects. We find a total Mi; m of 27.3 mag for the Gaus- 
sian and 27.2 mag for the Moffat profile (see FigureQ, where 
the total magnitudes are based on SExtractor's mag_auto 

6 The equation has a typographic error. The sum should be over i[ , 12 = 1 
to f'l , (2 = 3. 

7 The Moffat profile 1 Moffat 1969) is a modified Lorentzian with a variable 
power-law index that takes into account the flux in the wings of the intensity 
profile which are not included in a Gaussian profile. 
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FIG. 3. — Detection efficiency of the u image based on simulations planting 
objects of different intrinsic profiles into the image. The solid line represents 
a gaussian profile with a FWHM of l."3 and the dashed line represents a 
Moffat profile with a half-width at half maximum of l."2 and a power law 
index of 3.4. The u values are total magnitudes based on SExtractor's 
mag_auto apertures which are Kron-like (Kron 1980) elliptical apertures 
corrected for possible contamination. We find a total u\\ m of 27.3 mag for 
the gaussian and 27.2 mag for the Moffat profile, and for isophotal apertures 
which are more appropriate for LBGs, find a K tim of 27.7 mag for the Gaus- 
sian and 27.6 mag for the Moffat profile. 

apertures, which are Kron-like ( Kronj [1980|) elliptical aper- 
tures corrected for possible contamination. While total mag- 
nitudes are generally reported, isophotal apertures are more 
appropriate for LBGs, and yield a u\\ m of 27.7 mag for the 
Gaussian and 27.6 mag for the Moffat profile. This is one of 
the deepest M-band images ever obtained, with our sensitiv- 
ity being similar to those re ported in the Keck Deep Fields 
dSawicki & Thompsonl2005l) . 

This detection method does not match the detection method 
used in §3.4, and therefore does not constrain our detection 
efficiency of LBGs. As explained below, we use our prior 
knowledge of the positions of the sources which yields a dif- 
ferent completeness limit. However, this result gives the depth 
of our M-band image for comparison to other studies. 

3.4. Photometry through Template Fitting 

In order to obtain robust colors across images with varied 
PSFs, it is necessary to match apertures and correct for PSF 
differences. If the difference is minor, then methods to apply 
aperture corrections to account for the vari ations are appro - 
priate, such as the ColorPro software by ICoe et al. (2006). 
However, such algorithms don't perform well when the differ- 
ence in the PSF FWHM is large, such as when combining the 
high-resolution data from the HST (0."09 FWHM) with the 
low-resolution images obtained in this study with Keck (l."3 
FWHM). In this case, the uncertainties in aperture correc- 
tions are unreasonably large, and the low-resolution images 
are crowded such that objects overlap, making object defini- 
tions that are valid in both high and low-resolution images 
difficult to determine. 

In order to avoid these uncertainties, we use the TFIT 
dLaidler et al.l 120071) template-fitting method that uses prior 
knowledge of the existence, locations, and morphologies of 
sources in the deeper high-resolution UDF images to improve 
the photometric measurements in our low-resolution M-band 
image. This method creates a template of every object by 
convolving each object in the high-resolution image with the 
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PSF of the low-resolution image. These templates are then 
fit to the low-resolution image to determine the flux of all 
the objects in the M-band, relative to the flux in the UDF 
V-band image. We chose the V-band as the high-resolution 
reference image because it is closest in wavelength to the 
M-band, without being affected by the Ly-a forest over the 
redshift interval 2.5 < z < 3.5. The result is a very robust 
color which relates every object in the high-resolution V-band 
image to the low-resolution M-band, avoiding the problem of 
aperture matching between the two images while intrinsically 
correcting for the PSF difference. Using the V-band flux, 
the color is converted to a M-band flux, which inherits the 
same isophotal aperture as the high-resolution V-band image. 
The aperture correction used to obtain total fluxes for the V- 
band image is then also valid to obtain total fluxes for the 
M-band . For a more in-depth explanation, see lLaidler et al.l 
( 2007b. This technique is similar to others in the literature 
dFernandez-Soto et alJ1999tlLabbe et ai]|2005t IShaplev et al.1 
2005; Grazian et al. 2006), with the original version based on 
Papovich et all (2001 . 2001). 



We chose to use TFIT as it is publicly available, well doc- 
umented, and its performance is carefully tested. Like all the 
other methods, there are some constraints that had to be met 
to use this algorithm. The first is that the pixel scale of the 
M-band image must be an integer multiple of the pixel scale of 
the UDF V-band image. This was accomplished by drizzling 
the images such that the pixel scale of the M-band is 4 times 
larger than the V-band, as mentioned in §3.1. The second 
and third requirements are that the images must not be rotated 
with respect to each other, and the corner of the V-band image 
must coincide with the corner of the M-band image. These two 
requirements were met by rotating and trimming the V-band 
image using IDL. The resultant image was compared to the 
original image, and the difference in photometry was negligi- 
ble compared to the intrinsic uncertainties. We also improved 
our fit by source weighting the rms map before providin g it to 
the TFIT pipeline, as suggested by lLaidler et all (120071) . 

In order to avoid proliferating different catalogs with minor 
differen ces in object defi nitions, we adopt the object defini- 
tions of|Coe et al. (2006). These de finitions include t he cat - 
alogs of iBeckwith et al.1 (120061) and iThompson et all (2006), 
as well as detec tions performed on a white light image by 
( 20061) . By using identical object definitions as 
we can use the careful photometry for the 
B, V, i', z', J, and H bands already determined, and compare 
our redshift determinations knowing we have used identical 
apertures. 

TFIT requires a SExtractor catalog of the V-b and image 
as an input to the pipeline. The program sex s eg (Coe et all 
2006) was run on the segmentation map o flCoe et al.l (120 06 to 
provide the necessary inf ormation to TFIT while using the de- 
sired object definitions o flCoe et"ai1(l2006l) . The segmentation 
map defines which pixels belong to each identified V-band 
object. The sexseg program forces SExtractor to run us- 
ing a predefined segmentation map (for details, see ICoe et al.1 
2006). TFIT also requires a representative 2-D model of the 
M-band PSF, and is sensitive to the quality of its construction. 
We use the same Moffat profile fit described in §3.3 to model 
the PSF of the M-band image, and then use this as the transfer 
kernel by TFIT to convolve the V-band galaxy cutouts. We 
note that there is no significant spatial variation of the PSF 
across the field. 

In general, the higher the resolution and sensitivity of the 
high-resolution image, the better TFIT can model the sources 
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for the low-resolution image. If an input catalog is not com- 
plete enough, then unmodeled objects can act as an unsub- 
tracted background , slightly increasing the flux of all objects 
( lLaidler et al .1120071) . However, this has a limit, and eventually 
there are so many sources that are too faint to detect in the 
low-resolution image that galaxies are not well constrained 
given the substantial number of priors. This yields a large 
number of galaxies with unconstrained fluxes that increase the 
uncertainties of the nearby objects without yielding any new 
information. The UDF V-band image has substantially higher 
resolution and is deeper than the M-band image, and therefore 
a limit was put on the faintest galaxy used as a prior in TFIT. 
Only galaxies brighter than V = 29 mag are included in the 
input catalog to TFIT. The galaxies fainter than V = 29 mag 
are too faint to be constrained by the u image, and only add 
noise to the TFIT results. We stress that this is a conservative 
cut, and does not introduce an unsubtracted background. 

The quality of the resulting photometric fits can be evalu- 
ated through Figure|U which depicts four panels: the V-band 
image from the UDF, the u image from Keck, the model im- 
age, and the residual image. The model and residual images 
are diagnostics produced by TFIT, and are not used in the 
fitting process. The model image is a collage of the V-band 
galaxies convolved with the PSF of the M-band image, scaled 
by the TFIT flux measurement for each object. The residual 
image is the difference of the model and the M-band images. 
Ideally the residual image would be zero, but this is not the 
case (especially for bright objects), with multiple effects con- 
tributing to the imperfect residual. For instance, if the object 
in the V-band image is saturated, then it has the wrong profile 
for the M-band and leaves a residual. Alternatively, imper- 
fections in the modeled PSFs when scaled to large flux mea- 
surements of bright objects will also leave a residual. This 
effect was minimized by using a source-weighted rms map, 
although photometry of the brightest objects are imperfect. 
In practice, it is very difficult to perfectly align two images, 
the distortion correction is not perfect, and images generally 
have spatially varying PSFs. To minimize these affects, TFIT 
does a "registration dance", where it cross-correlates each re- 
gion of the model with the region of the data to find any local 
shifts. This registration dance was performed, which slightly 
improved the residual image, and leads to more robust pho- 
tometry. 

3.5. Sample Selection 

Our aim is to identify a sample of high-redshift galaxies 
that are suitable for constraining the star formation efficiency 
of gas at z ~ 3. A large fraction of objects for which fluxes 
are measured with TFIT are too faint to yield sufficient infor- 
mation regarding the object's redshift, and we therefore limit 
our sample to those objects with high S/N. We select objects 
based on their V-band magnitudes, since cuts in M-band would 
preferentially remove LBGs. The median M-band S/N of all 
objects decreases as a function of the V-band magnitude and 
drops below 3a at V > 27.6 mag. We adopt this V-band mag- 
nitude cut to include the majority of high S/N u objects, while 
removing S/N < 3 objects. We note that this is a conserva- 
tive 3cr cut since most LBGs won't be detected in the M-band 
reducing the overall median S/N. 

In addition to removing low S/N objects, we wish to re- 
move objects with photometry affected by nearby neighbors. 
TFIT can identify such objects with the covariance index di - 
agnostic that uses the covariance matrix (lLaidler et al.l 2007). 
During the fitting of an object's photometry as described in 
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FIG. 4. — Quality of TFIT photometric fits of the UDF V-band and LRIS n-band imaging. The images cover the overlap region of the central 1 1.56 arcmin 2 of 
the UDF, with North pointing up and East to the left. Clockwise from upper left: (a) the V-band image from the UDF; (b) the w-band image from Keck; (c) the 
model image; (d) the residual image. The model image is a collage of the V-band galaxies convolved by the PSF of the n-band image, scaled by the TFIT flux 
measurement for each object. The residual image is the difference of the model and the w-band image. Ideally the residual image would be zero, but this is not 
the case for bright objects (see §3.4). 



§3.4, TFIT uses the singular value decomposition routine to 
perform a chi-square (\ ) minimization. This yields a covari- 
ance matrix which is used to calculate uncertainties via the 
square root of the variance (the diagonal element), as well as 
the covariance (the off diagonal elements) of all objects in the 
fit. The covariance index is the absolute valu e of the ratio 
of the off-diagonal and the diagonal elements (Laidl eret alj 
120071) . The maximum value of the covariance index is saved, 
along with the corresponding object ID, and yields informa- 
tion about how an object's photometry is affected by its most 
influential neighbor. Objects for which this ratio is much less 
than 1 are generally isolated objects, while objects with large 
covariance index values can have unreliable photometry. Mul- 
tiple cuts are implemented to remove objects whose photom- 
etry have been significantly affected. First, all objects with 
a covariance index greater than 1 are cut because their mea- 
surements are not considered reliable. All remaining objects 
are kept if one of two conditions apply: either they have a 
covariance index less than 0.5, or they have a V-band flux 
greater than twice that of the nearest neighbor. This approach 
balances the desire for a large sample with the need to ob- 



tain reliable photometry. Lastly, we only consider objects that 
are detected in all four ACS bands (B, V, i', and z') to facil- 
itate and improve color selection in §4.2. This requirement 
removes most galaxies at z > 4 as they have little flux in the 
B or redder bands. 

Table[T]lists the 1457 galaxies that are left after the V-band 
magnitude, covariance index, and V-band flux ratio cuts. Each 
entry i ncludes the object ID (matching those from lCoe et alJ 
(2006)), R.A., decl.,' u magnitudes and uncertainties, and u 
S/N. It also lists information regarding the most significant 
neighbor, namely, its object ID, covariance index, V-band 
flux ratio, and separation distance. The u magnitude uncer- 
tainties include the uncertainties due to the V-band aperture 
correction made in ICoe et all d2006l) . as the accuracy of our 
u magnitudes depend on the accuracy of the V-band magni- 
tudes. Objects that are observed with a M-band flux less than 
3(7 significance are considered undetected in the M-band and 
are assigned a 3er upper limit. The magnitude limit is set to: 

m 3a =-2.51og 10 (3tT T FiT) + Z LRIS , (2) 

where Zlris is the zero point magnitude. The error distribu- 
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tion from TFIT is normally distribut ed and the u pper limits 
from TFIT are robust, as evaluated in lLaidlert 

4. PHOTOMETRIC SELECTION OF Z ~ 3 GALAXIES 

Star-forming galaxies early in their history, such as LBGs, 
exhibit a clear break in their Spectral Energy Distribution 
(SED) at the 912 A Lyman limit (Lyman break), as well as 
multiple absorption lines shortward of 1216 A by the Ly- 
man series. Photons bluer than the Lyman limit are not ob- 
served because the interstellar gas intrinsic to the galaxy and 
components of the foreground gas along the line of sight 
of the galaxy are optically thick at A < 912 A. At a red- 
shift of > 2.5, these spectral features are redshifted to opti- 
cal wavelengths, with the Lyman break entering the w-band 
(^3500 A). The ability to observe the Lyman break optically 
allows large samples of high-redshift galaxies to be identi- 
fied based on multicolor photometry where LBGs are selected 
by a strong flux decrement shortward of the Lyman limit, 
and a continuum longward of rest frame Lyq (1215 A) (e.g., 
ISteidel & Hamiltonl 1 9921; feteidel et alJl995lll996allbT) . In ad- 
dition, photometric redshift determination algorithms can es- 
timate object redshifts using galaxy SED templates 8 , which 
include additional information other than the Lyman break, 
such as the slope of the SED, the Balmer break at 3646 A, 
and the more pronounced 4000 A break, due to the sudden on- 
set of stellar phot ospheric opacity by ionized metals and the 
Call HK doublet dHamiltonll 19851) . Although color selection 
and photometric redshifts utilize the same SEDs for selecting 
Z ~ 3 galaxies, they each have their strengths and weaknesses. 
We use a combination of color selection and photometric red- 
shifts to create our sample of LBGs (§4.2). We provide a de- 
scription of the photometric redshift process, color selection 
method, and a catalog of all objects and their photometric red- 
shifts below. 

4. 1 . Photometric Redshifts 

Photometric redshifts (hereafter, photo-z's) are a well 
known and robust procedure to dete rmine redshifts of galaxies 
when spectra are unavailable (e.g. [Koolll985l; lLanzetta et alJ 
ll998HBemtezl2000tlCoe et alJ2006t IHildebrandt et alJ2008l) . 
They have the advantage over color selection that they take 
into account all the colors available simultaneously in \ 2 fits 
to template SEDs and yield more precise redshift informa- 
tion with clear redshift confidence limits. The photo-z's also 
sample z ~ 3 galaxies in regions of color space that color se- 
lected samples avoid because of low-redshift galaxies, and 
therefore can provide a larger sample. However, photo-z un- 
certainties do not always include systematic errors caused by 
variations and evolution of galaxy SEDs compared to SED 
templates, and possible mismatches of SED templates (see 
§4.1.2). Such systematic problems and the lack of a large 
spectroscopic sample make it difficult to characterize the con- 
tamination fraction of photo-z selected z ~ 3 galaxies (see 
§4. 1 .3). Nonetheless, they provide the largest sample of LBGs 
for study. 

For each galaxy, photo-z codes produce a probability distri- 
bution function, P(z), representing the probability of a galaxy 
being at any specific redshift. However, the P(z) can have 
multiple peaks, especially at z ~ 3 where there is a degener- 
acy with galaxies at z ~ 0.2, which then translates into large 

8 We refer to these templates as SED templates throughout the paper to 
distinguish them from the galaxy templates discussed in §3. 



uncertainties for the photo-z (Bemtez 2000). The introduc- 
tion of the M-band data helps resolve the photo-z degeneracy 
and improve the photo-z fits for z ~ 3 galaxies as it targets 
the most dominant signature in their SED, the Lyman break. 
We present photo-z's for the entire sample of galaxies with 
M-band data from §3.5 in Table |2] but caution against using 
them blindly to select galaxies at z ~ 3. We recommend mak- 
ing cuts on the sample to select galaxies with good x m od an ^ 
ODDS (for a description of these parameters, see §4.1.1). 

4.1.1. Bayesian Photometric Redshifts 

There are many diff erent photo-z codes available, and 
iHildebrandt et al.l J2008) explain the benefits of the different 
algorithms. We chose to use the Bayesian photo -z's (BPZ) 
dBenftezl200"0llBenftez etaLll2004tlCoe~e t al. 2006) to be con- 
sistent wit h past photo-z's d e termined for the UDF wi thout u- 
band data (ICoe et al.ll2006h . IHildebrandt et al.l (120081) advise 
using the SED templates supplied with their respective codes, 
since user-supplied SED templates can cause problems. How- 
ever, a re-calibration of the SED template set improves 
the performance of the photo-z r edshifts, and we u se the 
re-calibra t ed SE D templates from iBenitez et al.l (120041) and 
ICoe et all (2006) that have been extensively tested with BPZ. 
These re-calibrated SEDs are based on the star-forming (Im), 
spiral (Scd, Sbc), an d elliptical (El) galaxy templates from 
IColeman et alJ (119801) . the star bursting galaxy templa t es with 
different reddening (SB2, SB3) from Kinney et al] (119961) . 
and the faint blue galaxy SEDs with ages o f 25 and 5 Myr and 
metall icities of Z=0.08 without dust fromlBruzual & Chariot! 
(2003), described in section §4.1 of ICoe et al.l (120061) . We 
interpolate between adjacent galaxy SED templates for an 
additional two SED te mplates in t he ph oto-z fit, similar to 
IBenitez et al.1 42004) and lCoe et al.l (120061) . 

SED templates are not always a good match for each spe- 
cific galaxy, and when it is not possible to get a good fit to the 
SED template, then the resulting redshift may not be accurate. 
As a diagnostic of the goodness of fit, the BPZ code provides 
a reduced chi square (xl) value. However, high yQ values 
do not always indicate an unreliable redshift. Bright galaxies 
with small photometric uncertainties will have larger xl val- 
ues than faint galaxies with larger photometric uncertainties 
for the same numerator in \ 2 (also known as the variance), ye t 
have more reliable redshifts (see Figure 21 in lCoe et al.l2 006). 
This problem occurs because the systematic uncertainties of 
the SED templates are not taken into account, making xl no 
longer represent the relative quality of the fit. Nonetheless, 
a mechanism to evaluate the quality of the fits is r equired to 
trim th e sample to reliable redshifts. To this end, ICoe et all 
(2006) introduce a modified reduced chi square (x m od) va l ue 
that assigns an uncertainty to the SED templates in addition 
to the uncertainty in the photom etry of the galaxy . For clarity, 
we reproduce the equation from lCoe et al.l (1200 6) here: 



Xmod 



E 



(fa-fraY 

9 9 
CT f +CT f 

Ja JT 



(3) 



where f a are the observed fluxes, oy a is the error in observed 
fluxes, and fr a are the model fluxes, normalized to the ob- 
served flux£S 1 j7/ r represent the model flux errors, which are 
set by ICoe et all ?2006l) to a fr =max a (/ rQ )/15. While the 
definition of a/ T is arbitrary, it was picked such that the re- 
sultant x m od i s a more realistic measure of the goodness of 
fit. This is especially important for bright galaxies, as uncer- 
tainties in the templates dominate the error budget, and the 
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TABLE 1 
Catalog of u-band Objects 



ID a 


R.A. 


Decl. 


u (mag) b 


S/N 


Cov. ID C 


Cov. Index d 


Flux Ratio c 


Distance (arcsec/ 


1 


3 32 39.723 


-27 49 42.53 


24.69±0.01 


112.53 


12 


0.039 


0.007 


2.546 


7 


3 32 39.451 


-27 49 42.95 


25.76±0.02 


49.56 


10 


0.413 


0.112 


1.138 


8 


3 32 39.540 


-27 49 28.35 


26.18±0.04 


29.43 


80 


0.003 


0.003 


4.335 


13 


3 32 39.326 


-27 49 39.06 


26.71±0.05 


23.76 


16 


0.095 


0.060 


2.099 


14 


3 32 38.467 


-27 49 31.84 


25.28±0.02 


73.11 


56 


0.010 


0.012 


3.319 


15 


3 32 38.846 


-27 49 39.29 


28.92±0.34 


0.34 


22 


0.159 


0.591 


1.800 


22 


3 32 38.957 


-27 49 38.33 


28.91±0.35 


2.25 


19 


0.281 


12.129 


1.440 


24 


3 32 39.053 


-27 49 38.76 


27.34±0.08 


13.27 


22 


0.296 


0.911 


1.484 


33 


3 32 39.179 


-27 49 36.13 


27.07±0.06 


18.38 


48 


0.018 


0.393 


2.902 


35 


3 32 38.763 


-27 49 36.80 


27.73±0.11 


9.78 


19 


0.026 


10.179 


2.880 



NOTE. — Table |T|is published in its entirety in a machine-readable form in the online version of the Astrophysical Journal. A portion is shown here 
for guidance regarding its form and content. 

:i ID numbers from lCoe etal] 120061) . 

k Total w-band magnitudes for the same aperture as the V-band in Coe et al. i 2006). 

c ID of the most significant neighbor affecting the photometry according to the maximum covariance index. 

The maximum covariance index is the absolute value of the ratio of the variance and the covariance. 
c The ratio of object's V-band flux with the V-band flux of its most significant neighbor. 
^ Distance to the most significant neighbor. 



TABLE 2 

Catalog of Bayesian Photometric Redshifts 



ID" 


, b 
Zb 




0DDS d 




•^mod 


Zil 8 


f*l c 


0DDSl h 


1 


I7+ 011 


3.00 


1.000 


16.40 


0.20 


17* - 07 


3.00 


1.000 


7 


001 +o.io 
u - ul -o.oi 


6.67 


0.990 


7.59 


0.49 


Om+o.io 


6.67 


0.990 


8 


"'-"-0.15 


1.00 


1.000 


21.81 


0.01 


52 +<) 07 


1.00 


1.000 


13 


43+0.14 


3.33 


1.000 


8.00 


0.21 


43+0.11 


3.33 


1.000 


14 




2.33 


1.000 


46.40 


0.06 




2.33 


1.000 


15 

22 


52 +o.l5 
"■■"-a 16 

n T 3+0.41 
3 ""-0.41 


1.67 
5.00 


0.970 
0.999 


1.98 
7.34 


0.09 
0.44 


0.52t?i 

3 T 3+0. 19 
■'•'"-0.29 


1.67 
5.00 


0.520 
0.977 


24 

33 


1 3Q+0.78 
1,J *-0.23 

0.78^}? 


7.00 
3.67 


0.363 
1.000 


1.56 
1.96 


1.68 
0.08 


1 3Q+0.12 
n 7O+0.09 


7.00 
3.67 


0.247 
1.000 


35 




6.00 


0.996 


4.81 


0.60 


1 69+° 16 


6.00 


0.996 



zt>2 g 



t b 2 c 0DDS2 h 



i o i +0.02 
1.81_ a02 



6.00 1.000 



■r-7+0. 10 

J/ -0.02 


2.00 


0.443 


yu -0.10 


3.33 


0.023 


65 +0.22 
OJ -0.14 


6.67 


0.429 



02+ 03 
u ' uz -0.01 



6.67 0.004 



NOTE. — Table [2] is published in its entirety in a machine-readable form in the online version of the Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. In addition to the most likely redshift, we report the two most likely redshifts for each 
galaxy when available, along with the redshift ranges for each peak and the fractions of P(z) contained in those peaks. 
a ID numbers from lCoe et afN2006l) . 

Bayesian photometric redshift (BPZ) and uncertainty from 95% confidence interval. 
c Template SEDs used in the BPZ code as described in §4.1.1, where l=El_cww, 2=Scd_cww, 3=Sbc_cww, 4=Im_cww, 5=SB3_kin, 
6=SB2_kin , 7=25Myr, and 8=5Myr. Non-interger values are for templates interpolated between adjacent templates. 
d Integrated P(z) contained within 0. 1(1 +Zb)- 
e Chi square quality of photo-z fit. 

1 Modified reduced chi-square fit, where the templates are given uncertainties. 
g Redshift ranges for the two most likely peaks. 



Integrated P(z) contained within the local minima of P(z) for each peak. 

X\, values are not useful. The reported x 2 mo d values are re- 
duced chi square values, obtained by dividing by the number 
of degrees of freedom, v. The number of degrees of freedom 
is the difference between the number of filters observed and 
the number of parameters (in this case there are three fit pa- 
rameters, redshift, template, and amplitude). The minimum 
number of filters used is 5 and the maximum is 7, so the range 
of v in our study is 2 < v < 4. We note that Xmod * s calculated 
after the photo-z determinations, and does not affect P(z). 

If the quality of the fit to the SED template is good, 
then the ODDS parameter is useful in measuring the spread 
in P(z)- A galaxy with high ODDS has a single peak in 
P(z), while multiple or very wide peaks yield low ODDS. 
In general, restricting the photometric sample to those ob- 



jects with ODDS > 0.9- -0.99 yield clearly defined redshifts 
dBemtezlteOOOt iBenftez et~alll2004t ICoe et alj|2006h . In this 
paper we are conservative and restrict our sample to objects 
with the best vales of ODDS, those with ODDS > 0.99. Ad- 
ditionally, selecting galaxies based on SED template type 
(tb) can useful when selecting a specific type of galaxy, 
where l=El_cww, 2=Scd_cww, 3=Sbc_cww, 4=Im_cww, 
5=SB3_kin, 6=SB2_kin , 7=25Myr, and 8=5Myr. For in- 
stance, in selecting LBGs we constrain ourselves to galaxies 
with tb > 3, which only include star-forming galaxy templates. 
The redshifts, redshift uncertainties for a 95% confidence in- 
terval, t b , ODDS, xl, and x m0 d are a11 tabulated in Tabled 
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4.1.2. Photometric Redshift Measurement Uncertainties 

It is important to understand the origins of photo-z uncer- 
tainties in order to have confidence in their values. There 
are two types of uncertainties in photo-z's: 1) photometric 
meas urement uncertainties and 2) template mismatch vari- 
ance (lLanzetta et al.lfl998t iFernandez-Soto et al.ll2001l 120021: 
IChen et ail 12003b . Photometric measurement uncertainties 
are well understood, and they are responsible for the width 
of P(z), which determines the reported photo-z uncertainties. 
Faint galaxies with larger photometric measurement uncer- 
tainties will yield larger uncertainties in the photo-z's than 
brighter galaxies. If the photometric measurement uncertain- 
ties are very large, then the photo-z will be poorly constrained, 
because this results in multiple peaks in P(z) corresponding 
to many possible redshifts. The other possible uncertainty 
comes from template mismatches, which is a systematic er- 
ror due to the finite number of templates used in the photo-z 
determination. Not all galaxies will be well represented by 
our template SEDs, yielding large xt values. 

One method to decrease template mismatch errors is to in- 
troduce more template SEDs, since that increases the chance 
that there exist good matching SED templates for each galaxy. 
While this can improve low-redshift performance, it also in- 
creases the number of degenerate soluti ons and therefore 
gives poorer high-redshift performance ( Hild ehrandt et al.l 
2008). Degenerate photo-z's occur when different SED tem- 
plates fit the photometry equally well, resulting in multiple 
peaks in P(z) and therefore multiple possible redshifts. We 
consider galaxy redshifts degenerate if they have two or more 
peaks in P(z) at 95% confidence separated by Az > 1. In 
these cases, the resultant reported uncertainties are very large 
and the galaxy redshift is poorly constrained. We are mainly 
interested in good performance at high redshift and therefore 
do not increase the number of SED templates. 

Template mismatches are also the cause for "catastrophic" 
photo-z errors that occur, where the photo-z is incorrec t and 
the u n certainty does not include t he correct reds hift (lEllisi 
119971: IFernandez-Soto et al.l 119991: iBenftezI l2000h . Catas- 
trophic photo-z errors typically occur because multiple peaks 
in P(z) are incorrectly suppressed, leaving only one peak. The 
suppression of multiple peaks in bright galaxies likely oc- 
cur because of their small photometric uncertainties yielding 
large x\ values without taking into account the systematic un- 
certainties in the SED templates, which exaggerate the differ- 
ences between the different SED template fits that correspond 
to different peaks in P(z). This suppresses the peaks at other 
redshifts, resulting in possible "catastrophic" photo-z errors 
(Dan Coe, private communication). 

It is likely that the incorrect suppression of peaks in P(z) 
can be fixed through the introduction of SED uncertainties in 
the initial \ 2 fit, although such an addition requires an under- 
standing of those uncertainties using large surveys with both 
photo-z and spec-z's. This is neither in the scope of the UDF 
or this paper, and is being investigated elsewhere (Coe et al. 
in prep). In the mean time, the best we can do is reject photo- 
z's based on their Xmod vames - However, x mod * s not a * rue 
statistical test like x\n and therefore cuts normally appropri- 
ate for x\ are not valid for Xmod- Additionally, galaxies with 
J and //-bands data have median x mo d l al E er by ~ 0.7 than 
those without infrared (IR) data, although the inclusion of the 
IR data improves the reliability of th e photo-z's beca use of 
the significantly increased lever arm (ICoeetalJ 120061). We 
therefore don't use the same cut as used in ICoe et al. N2006T) 



(Xmod < 1)' but rather use Xmod to conservatively remove pos- 
sible bad photo-z fits. While we do present all the photo-z fits, 
we only include those with x mo d < 4 f° r our z ~ 3 galaxy sam- 
ple. This cut only reduces the total number of galaxies with 
photo-z's by ~5% (to 1385 galaxies). 

4.1.3. Comparison of Photometric and Spectroscopic Redshifts 

In order to test the accuracy of the photo-z's, we compare 
the redshifts with spectroscopic redshifts (spec-z's). We com- 
pile a list of 100 reliable spec-z's in the UDF that match our 
sample from §3.5 (see Table O 9 . In this sample, 18 spec-z's 
are fr om the VIMOS VLT Deep Survey (VVDS: lFevre et al l 
2004), where we only include redshifts with 95% confidence 
and multiple lines. Another 22 reds hifts come from th e 
GOODS VLT VIMOS survey (VIMOS : lPopesso et al.ll2009l) . 
where we include redshifts with A or B quality spectra. These 
spectra have good cross-correlation coefficients of the spectra 
with the templates and multip le lines are well iden tified. An 
additional 6 redshifts are from Szokolv et al. (2004) using the 
VLT FORS1/FORS2 spectrographs, where we use only those 
flagged as 'reliable' redshifts (quality flags "2" or "2+"). The 
remaining 57 redshifts are from the GOODS VLT FOR S2 sur- 
vey (FORS2: IVanzella et al.ll200l 120061 [20081 120091) . where 
we only include redshifts from A or B quality spectra. We 
do not include redshifts from the slitless spectra obtained as 
part of the Grism ACS Progra m for Extragalactic Science 
(GRAPES) dPirzkal et alJ 120041) . since the redshift determi- 
nations do not provide an independent check to photo-z's be- 
cause photo-z's were used to help identify the emission lines 
(IXu et alJl2007h . 

The photo-z's agree relatively well with the spec-z's (see 
left panel of Figure |5J, and have 100% agreement in the red- 
shift interval of interest (2.5 < z < 3.5), although only five ob- 
jects have spec-z's at these redshifts. There are clearly some 
galaxies that have incorrect photo-z's at lower redshift where 
the M-band does not sample the Lyman break. Of the 100 
galaxies from Table [3] 97 have Xmod < 4, of which 93 have 
ODDS > 0.99. The galaxy at a spec-z of 1.99 (ID 6834) has a 
Xmod ~ 17 and the resultant photo-z should be ignored. 

Only object 8585 meets our criteria, but has a significantly 
different photo-z than its spec-z and does not include the cor- 
rect redshift in its P (z). This object, with a spec-z of z = 0.3775 
dFevre et al.l 120041) and a photo-z of 1.45 ±0.24, is a case 
where the spec-z may be wrong. In general we tried to min- 
imize this possibility by selecting reliable spec-z's, but the 
spec-z's can still be wrong as seen in some comparisons of 
IFernandez-Soto et al.l d2001l) . In our case, the photometric 
redshift yields a good fit to the SED templates, and the fit 
to the SED template at the spec-z is a bad fit. The published 
spectrum shows that the spec-z is mainly determined by the 
Ha line. This could easily be confused with the Oil line for 
a galaxy with a redshift of z = 1.43, which would be consis- 
tent with our photo-z. This leaves one galaxy with a possible 
"catastrophic error", although it is not in the redshift interval 
of interest (2.5 < z < 3.5). 

4.1.4. Improvement of photo-z's with the Addition of the u-band 

A comparison of the redshift interval (2.5 < z < 3.5) in the 
two panels of Figure [5] shows a significant improvement in 

9 Most of these redshifts are based on observations made with ESO Tele- 
scopes at the La Silla or Paranal Observatories under programme ID(s) 66. A- 
0270(A), 67.A-0418(A),171.A-3045, 170.A-0788, 074.A-0709, and 275.A- 
5060. 
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FIG. 5. — Differences in spectroscopic redshifts (spec-z's) and photometric redshifts (photo-z's) compared to the 100 spec-z's described in §4.1.3, where the 
photo-z uncertainties are for a 95% confidence interval. The left panel plots the difference in spec-z and photo-z with M-band, and the right panel plots the same 
but for the photo-z's without the M-band. The quality of the photometric redshift fit is described by the Xm 0( j ( see §4- 1.1), where objects with good fits, Xm 01 j < 4, 
are marked with triangles, and objects with Xm d > 4 are marked with open crosses. There is one point in the left panel at z ~2 that has a large X^od v ~ ^ 
should be ignored as it is not a good fit to any SED templates with the M-band. 



TABLE 3 

Catalog of Galaxies with Reliable Spectroscopic Redshifts in the UDF 



ID" 


Survey 


-spec 


z„ b 


2 c 
A mod 


0DDS d 


V (mag) 


u-V (mag) 


u—B (mag) 


V-z' (mag) 


3088 


FORS2 


0.13 


006 +o.io 

u - UD -0.06 


1.91 


1.00 


22.85±0.00 


1.99±0.01 


1.16±0.01 


0.48±0.00 


5670 


VVDS 


0.13 


09 +0 11 


0.03 


1.00 


21.23±0.00 


2.00±0.00 


1.12±0.00 


0.54±0.00 


1971 


VVDS 


0.15 


u - u3 -0.03 


0.07 


1.00 


20.46±0.00 


1.48±0.00 


0.82±0.00 


0.37±0.00 


2974 


VIMOS 


0.15 


1 5+0.11 

-0.11 


1.51 


1.00 


24.08±0.01 


1.97±0.03 


0.99±0.03 


0.49±0.01 


5620 


VVDS 


0.21 


0.21«;J| 


0.75 


1.00 


23.42±0.00 


0.88±0.01 


0.43±0.01 


0.01 ±0.00 


3822 


VIMOS 


0.21 


0.18!« : g 


0.59 


1.00 


19.14±0.00 


1.91 ±0.00 


0.87±0.00 


0.70±0.00 


1000 


FORS2 


0.21 


021+0.12 


1.08 


0.97 


23.39±0.00 


0.93±0.01 


0.46±0.01 


-0.01±0.01 


5606 


VVDS 


0.23 


i 6 +0.11 
-0.11 


0.05 


1.00 


21.14±0.00 


1.79±0.00 


0.82±0.00 


0.55±0.00 


5491 


VIMOS 


0.23 


17 +0 " 
u.i/_ 12 


1.94 


1.00 


22.24±0.00 


0.64±0.00 


0.30±0.00 


-0.15±0.00 


7847 


VVDS 


0.33 


31+0.13 

UJ -0.13 


0.02 


1.00 


22.00±0.00 


2.97±0.01 


1.36±0.01 


1.10±0.00 



NOTE. — Table[3lis published in its entirety in a machine-readable form in the online ver sion of the Astroph ysical Jou rnal. A portion 
is shown here for guidance regarding its form and content. Redshift surveys ar e VVDS iFevre et al. 2004), VIMOS i Po pesso et alj 
!2009l) , ISzokolv et alj 120041) , and FORS2 IVanzella et all2005ll20Ml2008ll2009[). V magnitudes are total AB magnitudes, and colors 
are isophotal colors. All photometry other than the M-band are from Coe et al. i 2006). Nondetections in u-band are given 3cr limiting 

magnitudes. 

a ID numbers from lCoe etal] 120061) . 

Bayesian Photometric Redshift (BPZ) and uncertainty from 95% confidence interval. 
c Modified reduced chi square fit, where the templates are given uncertainties. 



Integrated P(z) contained within 0.1(1 +Zb)- 

the photo-z's with M-band data. The left panel depicts photo- 
z's with M-band data while the right panel depicts photo-z's 
without M-band data. The M-band helps prevent catastrophic 
redshift errors that can occur because of a similarity in th e col- 
ors of low-redshift galaxies and high-redshift gala xies. (lEllisi 
Il997t iFernandez-Soto et aUll999UBemtezl l2000). Usually, 
this degeneracy causes P(z) to have multiple peaks repre- 
senting both pos sible redshifts in the absence of M-band data 
dCoe et al.l2 006). As discussed in §4. 1 .2, however, sometimes 
secondary peaks are absent yielding incorrect redshift uncer- 
tainties and possibly incorrect redshifts. For example, three 
galaxies (IDs 830, 4267, 5491) with xLd < 4 and DDDS> 0.99 
have catastrophic redshift errors without M-band photometry 
and gain accurate photo-z's based on their spec-z's after we 
include the M-band (e.g. see Figure |6). As expected, the u 
data clearly improve the photo-z's of galaxies at z ~ 3 when 
compared to the spec-z's. However, we are constrained to dis- 



cussing small number statistics, and the true contamination 
fraction is unknown. 

We would ideally like to compare a larger number of z ~ 3 
photo-z's to spec-z's to get a better understanding of the im- 
provement of the photo-z's with the addition of the M-band, al- 
though such spectroscopic data are lacking in the UDF. How- 
ever, since the M-band samples the Lyman break of z ~ 3 
galaxies, and we know that the photo-z's with M-band are more 
reliable than those without, we can investigate the changes in 
the photo-z's. We therefore compare all the photo-z's with 
Xmod < 4 of galaxies with and without the M-band data in 
Figure [7] which highlights two effects. The first is that P(z) 
changes markedly for 125 galaxies, mostly in the redshift in- 
terval 2 < z < 3, where the M-band probes the Lyman break. 
Of these, 102 change their photo-z's from z < 1 to z ^ 2-3, 
whereas 23 switch in the other direction. This change is a 
result of the code selecting different P(z) peaks as the most 
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FIG. 6. — Object 830: an example of a galaxy with z sp e C = 2.44 (Popesso et al. 2009) that has a catastrophic redshift error before the addition of the u-band. 
The left-hand side figures show the photo-z fit without M-band with z=0.24±0.15 and the right-hand side figures show the fit with w-band with z=2.72±0.37. 
The rectangles represent the assumed uncertainty of the SED template used for x^j, blue points are the H-band data, green are the ACS data, and red are the 
NICMOS data. Filled circles are for detected and observed data points, and squares are for unobserved data points. The addition of the M-band rules out the 
possibility of this being a low-redshift galaxy, and corrects the catastrophic redshift. 

91, to 365 galaxies that have a photo-z in this redshift inter- 
val either with or without the M-band (including the 23 that 
switched to low redshift). Of these 365 galaxies, 161 galaxies 
either had their photo-z changed or the degeneracy removed 
with the addition of the w-band. This shows that the addition 
of the M-band significantly changed the photo-z's of the z ~ 3 
galaxy sample by ~ 50%. 

4.2. Color Selection 

Color selection is an efficient means to select high-redshift 
galaxies, and extensive research has been carried out to deter- 
mine the best color criteria to minimize the int erloper frac 



probable redshift for galaxies having multiple peaks . 

The second effect shown in Figure [7] is the removal of the 
degeneracy of z ~ 3 and z ~ 0.2 photo-z's, due to the removal 
of one of the peaks in P(z). One hundred seventy five galax- 
ies went from being degenerate to non-degenerate, and are 
marked by the red crosses in Figure [7] Not all of the degen- 
erate photo-z's are removed; however, 5 1 galaxies are degen- 
erate between these redshifts as marked by the blue crosses in 
Figure [7] of which 17 were not degenerate before the addition 
of the w-band. The new degeneracies occur when the w-band 
best fit redshift is different compared to the best-fit to the other 
bands. The old degeneracies that are not removed occur when 
the M-band does not conclusively rule out another template, 
often because they are faint. Figure [7] also includes galaxies 
with Xmod < 4 in the photometric redshift fits selected to be 
at z ~ 3 using the color selection method described in §4.2 
(below), which is useful when comparing the two methods. 
Figure |8] shows an example of both effects, where the w-band 
changes the photo-z and removes a secondary peak in P(z) for 
a high-redshift galaxy. 

Out of 1384 galaxies with Xmod < 4, there are 274 galax- 
ies that have photo-z's in the interval 2.5 < z < 3.5 without 
M-band. The addition of the M-band increases this number by 



tion fro m low- redshift g alaxies or stars , e.g. |Steidel et al 



( 1996a b, 1999, 12001 : [Adelberger etaD <H)04); Coo ke et al 



(2005). The color selection criteria used in these studies 
are based on predicted colors of model star-forming galax- 
ies at high redshift, which are then confirmed with spec-z's, 
that result in known contamination frac tions between 3%-5% 
dSteidel et al l 120031; iReddv et al.ll2008l) . Such low contami- 
nation fractions are achieved by avoiding colors where low- 
redshift galaxies reside. While color selection techniques do 
not provide a complete sample of LBGs, they do an excellent 
job of selecting galaxies in a specific redshift range, as evi- 
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• Objects with xl„ < 4 

Color selected and x^, < 4 
^ Degenerate without u-band 
x Degenerate with u-band 




2 3 
Photo-z with u-band 

FIG. 7. — Photo-z's without u-band data versus photo-z's with u-band data 
for objects with Xmod < 4 in both cases. The red crosses designate galaxies 
where the photo-z's are degenerate without u-band, and the blue crosses are 
for degenerate photo-z's with u-band. The degeneracy in redshift occurs when 
galaxies have at least two significant peaks in the P(z), usually due to the de- 
generacy in the colors of a low-redshift galaxies and high-redshift galaxies. 
The orange dots are those objects that are selected using the color selection 
technique in §4.2. We note although the photo-z's with u-band may be in- 
correct, the addition of the u-band made a large difference in the photo-z's of 
z ~ 3 galaxies. 



denced by their contamination fractions. While the UDF data 
provide an extraordinary data set, they also use a different set 
of filters than used in previous color selection studies, mean- 
ing we must define new color criteria for LBG selection. We 
therefore develop and calibrate new color criteria for selecting 
z ~ 3 LBGs using the same methodology. Since our motiva- 
tion is to generate a sample of galaxies to put constraints on 
the star formation efficiency at high redshift, we choose our 
color selection criteria to best minimize possible low-redshift 
interlopers (see below). 

4.2.1. Color Selection Criteria 



Usi n g the same approach as ISteidel et all (Il996albl 1 19991 
120031) . lAdelberger et all (l2004lh and ICooke et al.l (120051) . we 
derive galaxy colors by evolving different galaxy SED tem- 
plates to high redshift convolved with the total throughput 
of the different filters shown in Figure [2] We include galaxy 
SED templates consistent with our photometric redshifts de- 
scribed in §4.1, with galaxy S EP t emplates fro m Kinne vet alj 
(IT996l) . IColeman et al.l (fl980h . and lBruzual & Chariot! (120031 
In addition we us e a 2. Gyr Elliptical Galaxy from the 
iBruzual & Chariot! (2003) synthesis code (E2G) since it is 
quite different t han th e elliptical galaxy SED template from 
Cole man et al.1 (1 19801) and represents possible low-redshift 
galaxies we wish to avoid. We apply the ^-correction for 
different redshifts and correct for the op acity from the inter- 
galactic medium by using estimates from Madau(1995). The 
resultant colors and redshifts of the SED template galaxies 
are used to determine the appropriate color criteria to maxi- 
mize the number of LBG candidates at z ~ 3 while minimiz- 
ing the contamination from objects at other redshifts. We test 
multiple color-color planes and find that for our set of filters, 
z ~ 3 LBGs can best be selected in a (u - V) versus (V - z') 
color-color plane. Figure [9] plots the expected colors of dif- 
ferent model galaxies at different redshifts in the (u-V) versus 
(V-Z 1 ) diagram. 

The region defining candidate z ~ 3 LBGs is indicated with 



the dashed black line in Figure [9] In order to avoid select- 
ing low-redshift galaxies, this selection region excludes SED 
template colors for z < 2.5. The deviations from SED tem- 
plates and photometric errors will cause intrinsic scatter in 
the color-color plane. We therefore leave ~ 0.2 mag between 
our color selection and the low-redshift elliptical galaxy SEDs 
that cause the largest contamination for galaxies at z ~ 3. In 
addition to the cut in the (u - V) versus (V - z') diagram, we 
also apply secondary color cuts using the (u-B) color to im- 
prove our color selections by removing potential interlopers. 
We also include a cut on V-band magnitude, where the bright 
end does not remove any LBG candidates and the faint end is 
the V-band magnitude determined in §3.5 to keep the S/N of 
the w-band > 3a. The following conservative constraints are 
used to select LBG candidates: 



(u-V) > 1.0, 
(K-fi)>0.8, 
(V-z')< 0.6, 
3(V-z')<(«-V)-1.2, 
23.5 < V <27.6. 

4.2.2. Reliability of Color Selection 



(4) 
(5) 
(6) 
(7) 
(8) 



In order to test our selection criteria, we compare our se- 
lection of the 100 galaxies with reliable spec-z's as described 
in §4.1.3 and Table [3] These spectra allow us to test the ef- 
ficacy of selecting targets via the (u-V) versus (V -z') color 
plane (Figure ITOb. No z < 2.5 galaxies with reliable spec- 
z's passed our color cut confirming that our cut effectively 
excludes low-redshift galaxies. We note that there are three 
Z > 2.5 objects that do not meet our color criteria in Figure 
[TOl The z = 3.68 object (triangle , object 865) is classified as 
a quasar bv lSzokolv et al.l d2004l) due to active galactic nuclei 
(AGNs) activity. The z = 4.77 object with (V-z') is a bright V- 
band dropout galaxy, whose V-band magnitude (27.26) is just 
bright enough to remain in the sample before the color cuts. 
Finally, the z = 3.80 is also excluded by our color cut due to a 
decrease in the V-band magnitude, reddening the color. 

There are a number of objects that could contaminate our 
sample of z ~ 3 LBGs because their spectra are unusual and 
thus do not match our model SEDs. While at brighter magni- 
tudes the color selection criteria include stars and AGNs, for 
our V > 23.5 sample the color selection criteria are mainly 
conta minated by low-redshift galaxies at z < 0.2 (Red dv et alj 
120081) . Additionally, about 1 /3 of the Distant Red Galax- 
ies (DRGs) fall within the c olor selection sample of Stei- 
del (Ivan Dokkum et all 120061) . which is similar to our cri- 
teria. DRGs are galaxies at z >2 that have faint UV lu- 
minosities, hav e previously und ergone their episode(s) of 
star form ation dFranx et a l. 2003), and have stellar masses > 
1O U M (Ivan Dokkum et al.ll2004ll2006l) . However, given the 
DRGs' estimated space d e nsity of (2.2 ±0.6 ) x 10" 4 Mpc" 3 
(Ivan Dokkum et all 120061) . iReddy & Steidel (120091) conser- 
vatively determine that the fractional contribution would be 
~2% for UV-faint sources such as our sample. 

As discussed earlier, to get a sense of the number of possi- 
ble interlopers, we would like to have a large number of spec- 
tra to determine the contamination fraction. However, we only 
have a relatively small number of spectra (100) as shown in 
Figure[l0]and Table[3] Instead, we com pare our color cuts to 
the (Un-G) versus (G-R) color cuts of lSteidel et ail d2003l) 
for redshifts 2.7 < z < 3.4, and assume that our sample has 
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FIG. 8. — Object 97: An example of a galaxy with z S p CC = 2.69 iPopesso et al. 2009) that changes photo-z's with the addition of the H-band. The left-hand side 
figures shows the photo-z fit without u-band with z=0.24+g and the right-hand side figures show the fit with u-band with z=2.90±0.38. The addition of the 
u-band rules out the possibility of this being a low-redshift galaxy, and therefore P(z) only has one significant peak. 

similar contamination. In Figure[9] we determine the redshift 
that corresponds to Steidel's color cut for a given SED tem- 
plate in (JJn — G) and (G-R). We then determine the (w — V) 



and (V - zf) colors that this SED template has at this redshift 
and compare those resulting colors to our color criteria by 
marking them with filled black circles in Figure [9] The filled 
black circles are outside or near the edge of our color selection 
criteria, which shows that we are more conservative in our cut 
for the redder SED templates and abou t equivalently c onser- 
vative for the bluer ones as the cut in lSteidel et al.l (120031) . We 
therefore i nfer that our contam ina tion fraction is com parable 
to those o flSteidel etail d2003l) and lReddv etaT] d2008l) . where 
the contamination fraction of z ~ 3 LBGs is ~ 3% for objects 
23.5 < V < 25.5 mag. 

5. SAMPLE OF Z ~ 3 GALAXIES 

Photometric redshifts and color selection are both good 
ways to select z ~ 3 galaxies. Photometric redshifts have the 
advantage of creating a larger sample since they can measure 
redshifts in regions of color space that color selected samples 
avoid because of low-redshift galaxies. They also use more 
information than color selection, including all colors simul- 
taneously to constrain the redshift. However, while the error 
rate of the photo-z's is not well defined, color selection is effi- 
cient, has a clearly defined contamination fraction, and allows 



direct comparisons to other studies. The completeness of our 
LBG selection is limited primarily on the M-band depth be- 
cause the ACS bands from HST are deeper. In order to charac- 
terize this completeness, we compare our color selected LBGs 
with other studies in §5.3. To justify such comparisons, we 
compare the redshift distribution of this sample in §5.1, and 
investigate our uncertainties from cosmic variance in §5.2. 

In choosing our color selection criteria we chose to be con- 
servative and create a less complete catalog with a small con- 
tamination f r action similar to that of lSteidel et alJ (120031) and 
Reddy ^FaTI (120081) of 3%. Depending on the purpose, a 
higher contamination fraction is acceptable in exchange for a 
larger and more complete sample. The photo-z sample yields 
a more complete sample, and may even have a similar contam- 
ination fraction based on our spectroscopic sample, although 
due to our small numbers at the redshift of interest, it is not 
clear at this point. If the lowest contamination fraction pos- 
sible is needed, a subset of z ~ 3 LBGs that are both color 
selected and photo-z selected are the most robust candidates 
available. 

We present both samples of z ~ 3 LBGs in Table |H for a 
total of 407 candidates, along with their photo-z's and col- 
ors. We distinguish the samples by designating them as either 
color selected, photo-z selected, or both. The photo-z sam- 
ple consists of galaxies in the redshift interval 2.5 < z < 3.5 
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FIG. 9. — (u — V) versus (V — z') color-color diagram depicting different 
model galaxies at different redshifts, where the solid lines are z > 2.5, the 
dashed lines are z < 2.5, and the symbols represent redshifts of 2.5, 2.8, 3.0, 
and 3.2 which i ncrease with increasi ng (u—V) color. We include galaxy SED 
templates from Kinney et al. 1 1996) for star bursting galaxies with different 
reddening (SB2, SB3), from Coleman et al. (1980) for star-forming galaxies 
(Im), spiral galaxies (Scd, Sbc), and elliptical galaxies (El), In addition, we 
use and two faint blue galaxy SEDs with ages of 25 and 5 Myr and metallic- 
itie s of Z=0.08 (25My r 5My r), and a 2.0 Gyr Elliptical Galaxy (E2G) from 
the Bruzual & Chariot 12003) synthesis code. The area above and to the left 
of the dashed black line is the region selecting candidate LBGs. The filled 
black circles are the bord erline colors for eac h SED template that correspond 
to the z ~ 3 color cut by Steidel et al. 1 2003), as described in §4.2.2. 

that have % > 3, ODDS > 0.99, and x mo d < 4 ' and contains 
365 galaxies. Of the 42 galaxies not included in our photo-z 
sample, two have z ~ 2.2, 1 1 have z > 3.5, eight others have 
ODDS < 0.99, and 21 others have x^od > 4 - The color selected 
sample contains 260 galaxies, all of which have photo-z's with 
z > 2, with 258 that have z > 2.5 and 11 that have z > 3.5. 
However, the overlap of the two samples is only 216 galaxies. 
We show the final LBG selection in Figure QT| which plots 
all objects from Table Q] on a color-color diagram, with 287 
galaxies falling in the color selection region corresponding to 
constraints [4] [6] and|7]from §4.2.1. There are 27 galaxies that 
are in the selection area in this diagram, but are rejected by the 
(u-B) color (constraint|5j, leaving 260 galaxies that are color 
selected. The objects marked as blue stars that are selected by 
photo-z's but not by color selection are generally galaxies at 
z ~ 2.8 that are missed by our color selection criteria in order 
to avoid elliptical and low-redshift galaxies. 

5.1. Redshift Distribution 

In order to understand the redshift distribution of the color 
selected LBG sample, we look at the photo-z's that meet the 
color selection criteria and our photo-z criteria of x m od < 4. 
This leaves a sample of 235 galaxies that have a mean red- 
shift of 3.0 ±0.3 in the redshift interval 2.4 < z < 3.8 (see left 
panel in Figure [T2l . with a median uncertainty in the photo- 
z's of ±0.4. Additionally, we investigate the redshift distri- 
bution by adding the probability histograms P(z) of the in- 
dividual galaxies, which yields a similar result (see the right 
pa nel in Figure[l2|. Th e redshifts selected are similar to those 
of ISteidel et al.l (12003b . with ~80% of our sample matching 
their reported redshift interval 2.7 < z < 3.4. In fact, their dis- 
tribution is very similar to ours, with a number of their LBGs 
falling outside of this interval. Given the large uncertainties 
in our photo-z's, we conclude that our color selec ted redshift 
distribution is similar to that of Stei del etail ([2003) within our 
uncertainties. The similar redshift distribution of our color se- 



FlG. 10. — Comparison of our selection criteria to 100 reliable spec-z 
galaxies as described in §4.1.3. Object redshifts are binned and represented 
with different symbols, where crosses are for objects with z < 2.0, diamonds 
for 2.0 < z < 2.5, open circles for 2.5 < z < 3.0, triangles for 3.0 < z. < 3.5, 
and squares for z > 3.5. The dashed line depicts our color cut. 
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FIG. 1 1 . — All objects from table[T]on a color-color diagram, where the 
dashed line refers to the selection criteria for constraints [4] [6] and |7j from 
§4.2.1. There are 27 objects that are in the color selection area in this di- 
agram but are not color selected because of constraint \5\ from §4.2.1. The 
218 galaxies that are both color selected and photo-z selected are marked by 
red circles, the 42 galaxies that are only color selected are marked by green 
squares, and the 147 galaxies selected only by photo-z's are marked by blue 
stars, yielding a total of 407 z ~ 3 LBG candidates. 

lected sample justifies our comparison of the number densities 
of LBGs in §5.3. 

5.2. Cosmic Variance 

The UDF has a very small volume, with our overlap area 
consisting of 11.56 arcmin 2 , which in the redshift interval 
2.5 < z < 3.5 is a comoving volume of ~ 38000 Mpc 3 . A 
single pointing with a small solid angle and such a small vol- 
ume is likely to be affected by cosmic variance, yielding larger 
than Poisson uncertainties in the LBG number counts. It is 
important to estimate the cosmic variance effect in order to 
understand the systematic uncertainties for comparisons with 
number densities of LBGs in the literature. 

We calculate the cosmic variance using the code 
from [ Newman & Davis ( 2002) and the prescription from 
Adelberge r et alj d 2005l) to get the fractional error per count 
<j/N for our given volume in the redshift interval 2.5 < z < 3 .5 
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TABLE 4 

Catalog of Candidate z~3 LBGs in the UDF 



ID a 




v 2 c 

^mod 


0DDS d 


Type e 


V (mag) 


u-V (mag) 


u—B (mag) 


V-z' (mag) 


22 


T TO+0.41 


0.44 


1.00 


2 


27.38±0.03 


1.53±0.36 


0.96±0.36 


0.49±0.04 


76 


2 62 +fl36 


0.52 


1.00 


2 


26.47±0.02 


1.39±0.16 


1.19±0.16 


0.11 ±0.03 


84 


i i i +0.40 

■'■''-o^o 


0.01 


1.00 


1 


26.56±0.02 


2.40±0.35 


1.69±0.35 


0.03±0.04 


97 


T O7+0.38 


0.13 


1.00 


1 


25.02±0.01 


2.01 ±0.06 


1.48±0.07 


0.08±0.01 


99 


9 o 1+0.37 


0.32 


1.00 


2 


24.74±0.01 


1.96±0.07 


1.51±0.08 


0.52±0.01 


101 


2 66+° 36 
z - oo -0.36 


0.17 


1.00 


2 


27.38±0.02 


1.45±0.32 


1.09±0.32 


0.09±0.05 


131 


3 06+ 40 

-'• uo -0.40 


0.21 


1.00 


1 


26.95±0.02 


2.02±0.35 


1.43±0.36 


0.03±0.05 


209 


2 67+ 36 


9.03 


1.00 


3 


26.92±0.02 


1.27±0.34 


1.26±0.34 


-0.31 ±0.05 


213 


3 14+0.41 


0.18 


1.00 


2 


26.91±0.02 


1.95±0.34 


1.19±0.35 


0.56±0.04 


230 


t 7n +0.46 


0.00 


1.00 


3 


25.68±0.01 


3.08±0.35 


1.29±0.36 


0.44±0.02 



NOTE. — Table [4l is published in its entirety in a machine-readable form in the online version of the Astrophysical Journal. 
A portion is shown here for guidance regarding its form and content. Photometric redshifts should only be used if they have low 
Xmod an ^ g 00 ^ ODDS values as described i n §4.1.1. V mag nitudes are total AB magnitudes, and colors are isophotal colors. All 
photometry other than the u-band ar e fromjCoe et al. ( 2006). Nondetections in u-band are given 3cr limiting magnitudes. 
a ID numbers from Coe et al. ( 2006). 

Bayesian photometric redshift (BPZ) and uncertainty from 95% confidence interval. 
c Modified reduced chi-square fit, where the templates are given uncertainties. 
d Integrated P(z) contained within 0.1(1 +Zb)- 

c Type of LBG selection, where l=both color and photo-;: selected, 2=only photo-? selected, and 3=only color selected. 



. The variance is determined from the integral of the linear 
regime of the cold dark matter (CDM) power spectrum (P(k)), 

4dm = / P(k)\W(k)\ 2 d 3 k , (9) 

where W(k) is the Fourier transform of our survey volume. 
Since we want the variance of galaxy counts rather than CDM 
fluctuations, we need to correct for the clustering bias (b) of 
LBGs to get their variance (erg), where a 2 ~ b <Jq DM - The 
galaxy bias for typical LBGs is then calculated from the ra- 
tio of galaxy to CDM fluctuations in spheres of comoving ra- 
dius 8 /i _1 Mpc, where a%(z) represe nts the CDM fluctuation s 
for our redshift, and b = os^/o^z) ( Adelberger et al.ll2005l) . 
The resulting variance depends on a fit to the LBG correlation 
function £ g (r) = (r/r„)~ 7 , where r a is the spatial correlation 
length and 7 the correlation index. The galaxy variance is 
then 

72(r fl /8^F 
8 ' s (3- 7 X4-7)(6-7)27 

59.3) from which we can then calculate 



dPeeblesll 19801 eq. 
the fractional error per LBG count. 

Empirical fits to the correlation function yield differing val- 
ues for r„ and 7 depending on the sample, redshift distribu- 
tion, luminosity range, and redshift, which in turn affect the 
and our fractional uncertainty (TAdelberge r et alj 
2005; I lildebrandt et al.l I2007L 120091: iKashikawa et al.1 120061: 



Lee et al.ll2 006; Ouchi et al. 2004, 2005; [Yoshida et aU2008l) . 



We avoid looking at samples covering a small are a of the 
sky such as the HDF (Giavalisco & Dickinson 200l|), as such 
studies are also p lagued by cosmic variance as shown in 
lOuchi etail (120051) . The values in the larger studies gener- 
ally vary between r a ~ 2.8-5.5 and 7 ~ 1.5-2.2, with an in- 
creasing r and 7 with increasing luminosity, i.e., the brighter 
LBGs are more strongly clustered. There is also some minor 
evolution with r edshift, where the highe r redshift galaxies are 
more clustered (Hildebrandt et al. 2009). These values result 



in erg lbgs °f ~ 0.56 — 1.1 which correspond to a fractional 
error per count of ~ 0.14-0.28. Our sample includes the 
fainter less clustered LBGs at z ~ 3, so we are on the less 



clustered side of this range. We therefore adopt a fractional 
error per count (a/N) of ~0.2 in the rest of this study, which 
suggests that we could detect a relatively large over-density or 
under-density in our small volume. In fact, there is evidence 
of an over-density of z ~ 3.7 galaxies in the Chandra Dee p 
Field South, of which the UDF is a part dKang & Imll2"009h . 
However, no clear over-density is indicated in the correlation 
length measure d from the GOO Ds survey covering the same 
area on the sky dLee et al.ll 2006) at slightly higher redshift. We 
use the above estimate of the cosmic variance for our small 
field of view to constrain our results of the number densities 
of LBGs in §5.3. 

5.3. LBG Number Counts 

The number counts of z ~ 3 LBGs per unit of magnitude 
indicate the completeness of the LBG selection, and can be 
compared to number counts from other studies. Such data 
are only available for color-selected samples, and we there- 
fore only use our color-selected sample in our comparison. 
The same comparison could be accomplished with the lu- 
minosity function, which we do not calculate because we 
only have one pointing and our comoving volume is small. 
In other words, we have a small number of LBGs that in 
conjunction with the uncertainty due to cosmic variance dis- 
cussed above, would not yield meaningful constraints on the 
luminosity function. Additionally, to calculate the luminos- 
ity functio n, Monte Ca r lo sim ulations need to be run as de- 
scribed in Red dy et al.l ((2008) that are computationally pro- 
hibitive given our complex analysis tec hnique to get reliabl e 
photometry with largely varying PSFs dLaidler et alj l2007h . 
We therefore compare our results with the number counts 
from other studies. We stress that we are comparing number 
counts that are not corrected for completeness, and therefore 
the counts will fall at faint magnitudes in each study due to 
sample incompleteness. 

Our ground-based LRIS images have much greater PSF 
FWHMs (~1".3) than the HST V-band image (-0.09"), thus 
affecting the observed number counts. The dominant effect 
is caused by the blending of neighboring objects that affects 
the color of objects and therefore their selection. In addi- 
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Photo-z Photo-z 

FIG. 12. — Histogram of 235 color selected photometric redshifts with Xmod < 4. The mean redshift is 3.0 zb 0.3, and the median uncertainty of the photo-z's is 
zt0.4. The left panel is a traditional histogram, and the right panel is the sum of the probability histograms P(z) of the individual galaxies. The redshift distribution 
selected are similar to those of Steidel et al. (2003). 



tion, at significantly lower resolution, isolated compact faint 
objects have part of their flux lost to the noise floor of the 
background. This causes the faintest objects to go undetected 
in the low-resolution images, even if they would have been 
detected in a similarly sensitive high-resolution image. We 
use the HST V-band image to determine our high-resolution- 
detection (HRD) number counts and correct for this resolu- 
tion effect. We convolve the HST V-band image with the PSF 
modeled from the LRIS V-band image (FWHM of -1") that 
was taken concurrently with the M-band data. This yields a 
low-resolution-detection (LRD) image from which a segmen- 
tation map is generated using SExtractor. 

The final photometry is measured with this new segmenta- 
tion map in all bands using sexseg as discussed in §3.4, and 
then the same color selection is used as discussed in §4.2.1. 
Figure [13] shows the number counts per half magnitude bin 
per square arcminute fo r both the HRD a nd LRD, along 
with number counts fromlReddv et al. (2008), th e Keck Deep 
Field (KDF) (ISawicki & Thompsonll2006h . and ISteidel et al.1 
(1999). For the HRD sample, we include only Poisson un- 
certainties for reference, while for the LRD we include both 
Poisson uncertainties and the fractional error per count of 0.2 
to take into account the cosmic variance (see §5.2). The fig- 
ure shows that we are photometrically complete to V ~ 27 
mag, after which we start losing LBGs due to the sensitivity 
of the M-band image. In order to compare to other studies, we 
convert their 7?-band magnitudes from Keck to V-band magni- 
tudes from HST with a ^-correction at z ~ 3 for the Im galaxy 
template described in §4.1, which results in an R-V color of 
~ -0.15. The LBG number counts in Figure[l3]of the LRD's 
agree at the brighter end to within our uncertainties, however, 
at the fainter end (V > 26.0), our results are larger than the 
uncorrected KDF results, which is the only survey to probe to 
equivalent depths. 

Based on the turnover in the LBG number counts in Fig- 
ure [13] the KDF and our LRD study appear to be complete 
to V > 26 and V ~ 27, respectively. The cumulative num- 
ber counts for V < 26 is 3.7 ±0.6 LBGs arcmin -2 in our LRD 
study and 4.3 ±0.2 LBGs arcmin" 2 for the KDF, where the un- 
certainties are Poisson. These results are consistent with each 
other, without invoking cosmic variance, and suggest that the 
difference in number counts for V > 26 is due to differing 



completeness limits. 

In order to understand this better, we also include the KDF 
number counts corrected for LBG completeness in Figure [T3l 
Th is completeness correction is different than the one applied 
in (ISawicki & Thompsonl2006h . as it previously included the 
volume correction simultaneously. The correction applied 
here assumes that all z ~ 3 LBGs are well-represented by the 
colors of a fiducial LBG w ith (1) fixed age of 1 00 Myr (2) a 
fixed redshift z = 3, and (3) ICalzetti et al.l (1 19941) dust extinc- 
tion for E(B-V) = 0.2. The incompleteness is then calculated 
by planting objects with these colors in the KDF images and 
determining how many are recovered, with the uncertainties 
estimated via bootstrap resampling (Marcin Sawicki, private 
co mmunication). This is not as careful a correction as applied 
in (ISawicki & Thompsonll2006l) . but serves to investigate the 
completeness differences in our studies with respect to LBG 
detection. Our LRD study number counts are consistent with 
the KDF LBG completeness corrected number counts down 
to V ^27. This helps reinforce that we are likely complete in 
detecting z ~ 3 LBGs to V ~ 27 magnitude, making our study 
the deepest to date. 

6. SUMMARY 

We use newly acquired ground-based M-band imaging with 
a depth of 30.7 mag arcsec" 2 (1<t„ sky fluctuations) and an 
isophotal limiting M-band magnitude of 27.6 mag to create 
a reliable sample of 407 z ~ 3 LBGs in the UDF. We use 
the template-fitting method TFIT dLaidler et all2007l) to mea- 
sure accurate photometry without the need for aperture cor- 
rections, and obtain robust colors across the largely varying 
PSFs of the UDF ACS images (0."09 FWHM) and the M-band 
image (l."3 FWHM). The results are as follows: 

1 . We calculate photometric redshifts for 14 5 7 galaxies us- 
ing the Ba yesian algorithm of iBenftezI d2000l) . iBenftez et al.l 
(2004), and lCoe et all d2006l) . of which 1384 are rel iable with 
Y 2 „ H < 4. We find that the previous photo-z's by ICoe et all 
(2006) do a good job of determining redshifts even without 
the M-band if their uncertainties are taken into account. How- 
ever, these uncertainties are often quite large at z ~ 3 due to 
the color degeneracy of z ~ 3 and z ~ 0.2 galaxies. 

2. The M-band significantly improves z ~ 3 photo-z deter- 
minations: out of 1384 galaxies, 175 galaxies no longer have 
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FIG. 13. — Number counts of z ~ 3 LBGs in 0.5 mag bins per square 
arcminute not corrected for completeness. Magnitudes are slightly offset 
to make the figure more clear. The full number counts based on the high- 
resolution-detections are marked HRD and are purple circles. The LRDs 
are marked LRD and are in red squares. The LRD are the number counts 
that would have been measured if the object detections would have been 
made using ground-based images with FWHM of 1".2. For the HRD sam- 
ple, we include only Poisson uncertainties, while for the LRD we include 
both Poisson uncertainties and the cosmic variance esti mate (see §5 . 2). Ou r 
number counts are plotted along with the studies o f Reddy et al.!| 120081) . 
ISawicki & Thompson 12006), and Steidel et al. 1 1999). The only points cor- 
rected for completeness are the gold crosses for the [Sawicki & T hompsor] 
12006) data. 

degenerate photo-z's, and 125 of the galaxies changed their 
primary photo-z with the addition of the M-band. In fact, the 
addition of the M-band changed the photo-z's of ~ 50% of the 
z ~ 3 galaxy sample. 

3. We find that even when using the M-band photometry and 
restricting the sample of photo-z's to those with good x^ od , 
catastrophic photo-z errors can still occur, although they are 
rare (only 1 out of 93 galaxies with spectroscopic redshifts 
and photometric redshifts with xLa < 4 and 0DDS > 0.99). 
We found no catastrophic photo-z's in the redshift interval of 
interest, (2.5 < z < 3.5), although only five objects have spec- 
z's in this interval. In contrast, three galaxies at 2.5 < z < 3.5 
had catastrophic photo-z's before the addition of the M-band. 
The contamination fraction of the z ~ 3 photo-z sample is 
likely small as we sample the Lyman break for these galaxies, 



and show excellent spectroscopic agreement. However, given 
the small numbers, the overall error rate of the photo-z's is not 
well defined. 

4. We find excellent agreement of our color selected sam- 
ple with the spectroscopic z ~ 3 sample, with no low-redshift 
galaxies falling in our color selection area, confirming our 
chosen criteria. We specifically ch oose a conservative color 
criteria that are similar to the cuts of Ste idel etail (120031) such 
that we can inf er that the ~ 3% co nta mination fract i on for 
z - 3 LBGs of lSteidel etaLI (120031) and lReddv et aH d2008l) 
applies to our data set. 

5. The completeness of our LBG selection depends largely 
on the M-band depth because the ACS bands from HST are 
deeper. In order to characterize this completeness, we com- 
pare our color selected LBGs with other studies and find that 
we present the deepest sample of z ~ 3 LBGs currently avail- 
able, likely complete to V ~ 27 mag. 

This reliable sample of z ~ 3 LBGs can be used to further 
the studies of LBGs and star formation efficiency of gas at 
z ~ 3 through the most sensitive high-resolution images ever 
taken; the Hubble Ultra Deep Field. 
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